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ABSTRACT 


Hydroacoustic communication has not been used ex- 
tensively at ranges in excess of ten miles. A study of 
t1e shallow-water and surface ducting channels indicate 
that reliable sixty word per minute communication at ranges 
of fifty miles is possible under a variety of conditions, 
With reasonable power levels and equipment complexity. lhe 
channel can be approximated by a Rayleigh fading signal with 
additive Gaussian noise. Five frequency multiplexed binary 
channels are compared with a fifteen channel system utilizing 
frequency diversity, and with a fifteen channel system util- 
iging coding. A method of calculating the diversity error 
for components with unequal signal and or noise levels is 
worked out. Frequencies must be above 900 cps if shipping 
noise is significant, and below 1500 cps to minimize at- 
tenuation. +0 optimize several transmitter~-receiver para- 
meters, such as the signal interval, more long range ex-~ 
perimental data is required. 
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INTRODUCTION 


Except for a few very short range or very low rate 
systems, the communications capability of ocean acoustics has 
not been exploited. Voice modulation has been used extensively 
for ranges of a few miles and CW at slightly lonser ranges. 
Single explosive charge signals have been sent over a thousand 
miles via deep ocean channels. Aside from sonar these few 
examples cover most of the communications development. With 
the increasing interest in nuclear submarines, oceanography, 
and the industrialization of underwater real estate, the use- 
fulness of a reliable system which does not depend on radio 


or cables is apparent. 


In this thesis a single such problem will be con- 
Sidered, that is, communication between two stations in rela- 
tively shallow water using transducers near the surface. A 
range of fifty miles and a data rate of sixty words per minute 
will be used to evaluate various systems. If such a system 
were feasible it would provide a good link between ships ope- 
rating in the same geographic area, or between an oceanorraphic 
vessel and its remote data collecting stations, or any number 


of other purposes. 


Part I of this thesis will be a simplified description 
of the shallow water acoustic channel, which will point out 
the varameters important to a communications system. These 
include ambient noise levels, transmission losses, and channel 


stability. ‘these parameters will be used to arrive ata 





random process model of the channel, 


In Part II this model will be used to develop a com- 
munication system. First a single binary channel will te 
developed. this will then be extended as required to obtain 
the desired data rate at a range of fifty miles with an ac- 
ceptable error rate. Y%ince there are many Ways to handle the 
Peeolem Only one or two representing current technology will 


be discussed. 


In Part III the feasibility of hydroacoustic communica- 
tions under various ocean conditions Will be evaluated. The 
assumptions used in the first two parts and tneir effects will 


also be evaluated at this point. 


Part IV will summarize the results and indicate areas 


requiring additional researcn and or development. 


rary 


SHALLOW -wATER HYDROACOUSTICS 


1. General. In shallow water there are two principle methods 
of acoustic transmission. JIne first type is the surface channel. 
In tnis channel a large portion of the acoustic energy is 
trapped in a duct near the surface. Transmission in this 

duct is by surface reflection and upward refraction. ‘the second 
type is transmission by both surface and bottom reflections. 

The particular type of transmission which will occur is pri- 
marily dependent on the vertical profile of the velocity of 
sound in the ocean area. ‘the velocity of sound in sea water 
varies primarily with temperature and to a lesser degree with 
pressure and salinity. Over long ranges it is reasonable that 


both types of transmission may occur for part of tne distance. 


2. Tne Surface Channel. “the density of water decreases with 
temperature. tnis means that in general the water near the 
surface is Warmer. However, if a cooling condition exists, 
the surface water may be cooler or at the same temperature 

as deeper Water. This inversion causes the surface channel to 
exist. The velocity of sound increases with increasing pres- 
sure and with increasing temperature (1). Therefore, an 
inversion or isothermal layer will have a positive velocity 
gradient, i.e., velocity of sound increases with depth. Ac- 
cording to Snell's iaw, acoustic waves in this layer will be 
bent upward, reflected from the surface and bent upward again. 


Figure 1 is a diagram of this type of channel. It has been 
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Figure 2. Spherical and cylindrical spreading losses. 





Shown and experimentally verified, that a minimum layer depth 
is required to trap a given frequency signal (2). tnis depth 
Ss 


Imin = 2.7 x 107 (fg)7173 


i 


frequency (cps) 
g = duct velocity gradient (sec.7}) 
For a 1000 cps signal, and a 0.020 meter/second/meter velocity 


gradient this minimum layer depth is 304 meters. 


There are several types of transmission losses which 
mast be considered in this channel. These are: spreading 
losses, reflection losses, and absorption losses. ‘the spread- 
ing loss will be broken into spherical and cylindrical com- 
ponents. tlhe energy leaves the source initially as a spherical 
Wave. Beyond a certain point the energy is essentially trapped 
in the layer and the spreading loss becomes cylindrical. For 
purposes of calculation the transition point will be taken as 
half the horizontal range of a ray leaving the source at the 
layer depth horizontally and just reaching the surface. ‘+his 


is shown in figure 2. 


The second tyve of loss to be considered is the surface 
reflection loss. ‘the reflection coefficient for the water- 
air interface is approximately minus one. ‘This is true for 
calm water, but waves cause a certain amount of scatter loss. 
This loss coefficient has been determined experimentally (3). 
Tie following formula describes the loss function for a fre- 


quency-waveheight product less than 3: 





&= -10 log, (1-0.149 (fH)3/2) ab per bounce 


a) 
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frequency in kc/sec. 


oe 
il 


Waveheight in meters 


The third tyoe of loss to be considered is absorption. 
This takes into account primarily the thermal and viscous 
forces Which are otherwise neglected in deriving the wave 
equation of acoustics in a liquid medium. For the ranges of 
interest this loss is negligible below frequencies of 2 kc. 
A suitable formula for absorption losses which is accurate 


from frequencies of 2-25 ke is (4) : 


“oe ; 
“a= 7200 x 10 Ef + Birmeley 10-8 f< nepers per Meter 
ASO ee 


this sives a transmission loss of: 


L =7.280 x 103x «, db/kn. 


If it is assumed that the velocity of sound does not 
vary in the channel except linearly with depth, some rather 
simple equations can be derived which specify the geometry 
of the path in depth and range. SBeferring to figure 3, 
these formulas are (5): 

a = transducer depth (meters) 

L = layer depth (meters) 

R = horizontal range (meters) 

a = velocity gradient (sec.71l) 

€= 90° - angle of incidence 

Y = vertical depth used for calculation (meters) 


Co= lowest velocity of sound in the layer (meters/sec.) 


Co velocity R range 





depth 


Figure 3. Surface channel ray path calculation geometry. 
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Figure 4, The SOFAR channel. 
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If we now let Y = L (transducer at the surface), and G= 0 
then 2R is thefmax imam distance between reflections, and the 
minimum number of reflections over a distance D is: 


minimum number of reflections = D_ 
ZR 


Since the water devth does not enter any calculations 
for the surface channel, it is also a suitable deep water model 


as long as it is uninterrupted. 


3. TDne Bottom Reflection Channel. The bottom reflection 
channel only exists in shallow water. In very deep water the 
temperature reaches a minimum at some intermediate depth and is 
isothermal down to the bottom. This gives rise to a minimum 
velocity of sound, which penerally occurs at depths of about 
two thousand feet. Sound waves are reflected upward below 

this depth, and downward above this depth. ‘his is the SOFAR 
channel shown in figure 4 It will not be considered in this 
thesis. For the shallow water channel, eee depths up to 

600 meters will be considered. This includes about ten per 


cent of the oceans, mostly along the continental shelf (6). 


in shallow water the SOFAR channel rarely exists, and 
the effects of both bottom and surface reflections are signi- 


ficant. Some typical linearized velocity profiles are shown 


‘\ 





in fieure 5. For uniform velocity gradients the ray path 
equations as described for, surface cnannel may be used for 
the bottom reflection channel path calculations. For a, b, 
and c of figure 5, these calculations would have to be made 


in several steps. 


Except for the addition of a bottom reflection loss 
the bottom reflection channel losses are the same as those of 
the surface channel. Because of the greater channel depth it 
is obvious that the svherical spreading loss will extend for a 
greater distance. ‘his will be computed as R vice R/2 ce 
surface channel. In general, experimental data shows that 
the total losses for this type of channel approximate spherical 
losses for the whole transmission distance. ‘he surface channel 
has total losses somewhat lower than spherical spreading. 
The bottom reflection loss is a complex function of the bottom 
material and its structure. For grazing angles less than 
fifteen degrees however, the losses are small. A good expe- 
rimental figure is 0.3 db per bounce (7). For larger grazing 
angles the loss becomes quite large. Figure 6 is a transmis- 
Sion path diagram of the bottom reflection channel. Note that 


rays which graze the bottom at about 15 degrees, leave the 


source at a much smaller angle. 


Sample calculations are made for both the surface 
channel and the bottom reflection channel in Appendix A. 
Several computations are plotted in figure 7 together with 


some experimental data for comparison. Im all cases the 
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Figure 5. Shallow water velocity profiles. 
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Figure 6. The bottom reflection channel. 
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Figure 7. Sound pressure level at 50 miles as a function 
of frequency. All data corrected to a source level 
- of 179 db. at a radius of one meter. References for 
experimental data in parentheses. 
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sources were corrected to a sound pressure level of 185 db 
(ref. 0.00002 newt./m¢) at a 1/4 meter radius. 


SPL = 20 log) g 2 
Lo 
This is low enough to permit simultaneous transmission over 


each of ten channels. 


4, Noise. Ambient sea noise is a complex array of processes 
which are not all fully understood (13). Most noise sources 
are very frequency dependent. Below 1000 cps ship and traffic 
noise is the most significant. I1t usually peaks between 300 
and 800 cps and falls off rapidly at higher frequencies. For 
all frequencies of interest above 1000 cps the ambient noise 
is primarily dependent on sea state for all but calm Weather. 
The primary sources of this noise are thought to be wind, spray, 
cavitation, and precipitation. Other sources are marine life, 
molecular and thermal agitation and earthquakes. Table l 
gives the average relations between wind force, wind velocity, 
Wave height, and sea state. Yigure 8 is a plot of the average 
ambient noise in the ocean versus sea state. The highest ex- 
pected values are about 3 db above the average. A possible 
explanation is an irregular gusting wind. It should also be 
noted that deep water noise levels are generally 5 db lower 


than shallow water levels. 


If a transducer is to be carried on a ship, the ship's 
self noise is very important (14). At low speeds this will 
correspond to the ambient noise at sea state 2, but it may 


contain some high level discrete frequencies. At speeds above 
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Figure 8. Ambient noise in the ocean(13). Specific 
shallow water data indicated as follows: 
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about fifteen knots, cavitation noise can become much greater 


than the ambient noise, 





Table 1 
Sea Beaufort Wind Speed Wave Height Sea 
State Wind knots feet Description 
force (meters/sec. ) (meters) 
2 3 7-10 1-2 Scattered whitecaps 
( 325- 5.0) ( o3- oy) 
3 4 11-16 BAS Frequent whitecaps 


( 5.5= Guo (o-) 6) 


4 5 17-21 5-10 Many whitecaps 
( 8.5=-11.0) (1.5-3.0) 


5 6 PPAF, S17 Whitecaps everywhere 
(11.5-14.0) (2.5-5.0) 


6 2 28-33 12-26 Heaped up sea, 
(14.5-17.0) (3.5-8.0) blown spray 

5. Source Power. The output power of a transducer is limited 
by cavitation and acoustic non-linearity. Cavitation occurs 
when the local pressure at or near the transducer is less 
than the vapor pressure of the water. A sample calculation 
in Appendix B shows that this limits the transducer toa 
sound pressure level of 195 db (ref. 0.00002 newt./m-) ata 


transducer depth of twenty feet. 


The source level of the transducer is limited by the 
approximations required to obtain the acoustic wave equations 
from the equations of fluid dynamics. ‘he governing ap- 
proximation comes from the following total derivative (17): 


du _du, udu 


— 


G2 Wee ar 
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u = fluid velocity vector 

tC = time 

xr = fluid particle position vector 
if the second term on the right can be neglected, then the 
Wave equation can be satisfied. This is worked out in Ap- 
pendix B for spherical waves. The result is that tne sound 
pressure level at one meter radius must be less than 254 db 


iret. 0.00002 newt./m-). 


6. Random Process Channel Model. An input-output relation 
for the hydroacoustic channel will now be developed. The 
parameters which will be considered are: signal energy, noise 
energy, frequency dispersion, time dispersion, and multipath 


mecepvtion. 


Since the ambient sea noise is due to a large number 
of small sources, it will be assumed to be a zero mean, 
Gaussian random process with spectral density Ne. Np is then 
the noise energy density. ‘his noise is indevendent of the 
message sent. Therefore the received signal will be a 


message filtered by the cnannel plus additive Gaussian noise. 


This channel filter function is unknown however. To 
arrive at a channel description, two assumptions will be 
made. First, along any given transmission path, the channel 
is a pure attenuator. Second, tne receiver receives the 
Signal over a large number of paths, With tne received pnase 


for each path being a uniformly distributed random variable. 


bs 
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Wozencraft and Jacobs show that this description leads to 
a random process channel model, with the following input- 


output relationships (18): 


s{t) = m(t) V2" cos Wot input 
r(t) = a(t)m(t) V2’ cos (wWot-6)+n(t) output 
in which 
p(a) = 2ae7a* eyes 0, 
p(o) = 1_ 0490427 
27r 


The amplitude a(t) is thus a Rayleigh-distributed random 
variable. Experimental data verifies this model in many 
instances, but other distributions have also been observed 
(19,20). For the work to follow, a(t) will be assumed to be 


Rayleich-distributed. 


Before a communications system can be developed, 
frequency and time dispersion must be discussed. Frequency 
dispersion will limit the bawi length. Experimental data 
shows that in general, for ranges of interest, frequency 
dispersion is less than 2 cps (21). It is expected that the 
system will use a baud corresponding to a much larger band- 


width than this. 


Time dispersion of a signal is the range of delays 
observed over paths of different length. It limits the rate 
at which signals may be sent over conventional systems. 
there is very little experimental data on shallow water time 


dispersbon. The calculations in Appendix C show that about 


ft 





a4 


1/5 second is the lonsest time that need to be considered, 
at a range of fifty miles. This is about 40 percent of the 
dispersion for the SOFAR channel which is severely limited 
in this respect (22), and it corresponds to an experimental 


value observed in the Florida Straits (23). 


a 
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COMHUNICATIONS SYSTEM 


1. Frequency Selection. The optimum frequency is the one 
Which allows the greatest bandwidth and the signal to noise 
ratio. Neglecting the noise, low frequencies are the best. 
Oopreading losses are about the same for all frequencies. 

For a state 3 sea,surface reflection losses are 0.5 db./ 
bounce at 800 cps, 0.95 db./bounce at 1000 cps, and 2.0 dab./ 
bounce at 1500 cps. This becomes worse for higher frequencies 
or higher seas. Attenuation also increases with frequency 
and is first significant for the fifty mile range at about 
2000 cps. Therefore, primarily because of surface losses, 

an upper limit of 1.5 ke will be used. A lower limit of 800 
cps will be used. At low frequencies it is possible to get 
strong harmonics from shipboard machinery. this is general- 
ly worst at about 500 cps. One prime example of low frequency 


noise is propeller harmonics which may be very strong. 


A plot of noise and fifty mile signal strength is 
shown in figure 9 as a function of frequency. This plot 
uses the same data as figures 7 and 8. A noise level cor- 
rected to a 10 cps bandwidth is also shown. ‘nis is suf- 
ficient for the correlation time of 1/6 second which will 
tentatively be used. Experiments show that at/long ranges, 
the signal will remain about 80% coherent by autocorrelation 


for up to 1/4 second (23). 
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<¢. Binary Sisnalline. A simple binary communications 





System will now be developed for the shallow water channel. 
The signal set will be an orthogonal pair of cosine waves, 
go Ol lows: 

S_tt) = V2" cos Cou 

si (t) = ~2 cos Wit 
This set is considered the easiest to use for this system. 
Signals which are time orthogonal are not well suited to this 
channel since the time of arrival is unknown. Since both 


phase and amplitude are random processes, neither phase nor 


amplitude modulation was considered suitable. 


When the signal s,(t) is sent, the received signal is 
Gmmcne form: 


N 
a) oan cos (Wott+0;) + n(t) 
= 


where the i's represent the different multipaths. As noted 

in Part I, this leads to the following random process model: 
r(t) = a(t) cos(Wot + 8) + n(t) 

in which a and 9 are random variables with the following 

probability density functions at any instant of time: 


2a @7a* 0a 


p(a) 
p(9) 


This model and the following assumptions will be used to 


ee 04027 


arrive at a suitable binary receiver: 

a. a(t) and O(t) are constant over the signal inter- 
val T, but vary randomly from interval to interval. 

b. n(t) is a white Gaussian noise random process, 


of spectral density No, and it is independent of the signal. 





roel 


c. The two signals, s,(t) and S,(t) have equal 
energy and are equiprobable. 


ad. There is no intersymbol interference. 


A maximum likelihood receiver will be used. For this 
receiver, that message is selected for which the probability 
of the received signal conditioned by the messages is highest. 
The details of this receiver have been worked out by Turin 
(26), and for more general time-varying cases by Kailath (27) 
and Price (28). In the case of the narrow band signals chosen 
for the signal set, tne receiver was determined to be a de- 
Vice which correlates, squares, and sums the quadrature con- 
ponents before making a decision. The correlation receiver 
is shown schematically in figure 10. Figure 11 is an envelope 
detector receiver which can be shown to be essentially the 


Same e 


The correlation receiver can be explained heuristic- 
ally by a geometric argument. Tne quadrature nolse con- 
ponents of the received waveform are zero mean Gaussian ran- 
dom variables of variance N,o/2, where No is the mean square 
noise power. The quadrature signal components of the received 
waveform are also zero mean Gaussian random variables, but 
with a variance of E,/2, Where Ep is the mean square signal 


power. If a signal is present the received quadrature com- 
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Figure 10. Binary correlation receiver. 










select 





largest 
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Figure 11. Binary envelope detector receiver. 
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ponents are therefore zero mean Gaussian random variables 
with variance ak. These Gaussian random variables are 
illustrated in figure 12, in the s,-plane and the S,; -plane 
for the case when s,(t) is sent. In each case the amplitude 
of the resultant vector is a Rayleigh distributed random 
variable. 

The error probabilities for the likelihood receiver are 
developed in Wozencraft and Jacobs (29), by computing it 


first for the random phase case only, to obtain: 


-/2N 
Ec) = +6 


f 


Averaging this error over a Rayleigh amplitude distribution, 
the following error probability is obtained for the random 


phase and random amplitude channel: 


PE) = say 
Using this formula, the efficiency of the simple binary 
channel may now be determined. Referring to figure 9, with 
a 10 cps bandwidth and the weak signal line, table 2 summarizes 


the binary channel error probability at various frequencies. 


lableme 
Frequency E/N (db. ) e/N ele) 
800 a 790 CG. 0019 
1000 25 200 0.0050 
1250 > SENS 02030 
1500 8 org) Olle 


3. Channel Capacity and Coding. In 1948, C.#. Shannon 
showed that whenever signals are transmitted at some rate, 
R, less than channel capacity, C, then an arbitrarily small 


error rate can be achieved by coding with a sufficiently 
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Figure 12. Received signal in the s,-plane and Ss, ~piane 
when message s,(t) is sent. 
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Figure 12A. The binary symmetric channel. 
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lons block length (30). For the binary symmetric channel 
With crossover probability bp, as shown in figurel2A, the 
channel capacity has the rather simple form (31): 

Ges VAT [ 1 + plogep + (1-p) logs (1-p)] bits per second 
where T is the signal time required to achieve p. It can 
also be shown that for n arbitrary channels being used in 
parallel, the capacity is: 

C2 £0, 

= 


+he 
where,Cj are the individual channel capacities. This last 


fFOrmula holds for simple parallel use in which a decision 

is made on each channel. If the decision is made on the col-~ 
lective output probabilities of the individual channels, then 
the channel capacity is somewhat greater, since in the first 


method some information is lost by the intermediate decision. 


A sixty word per minute data rate requires approxima- 
tely thirty bits of information per second. If an average 
channel capacity of five bits per second can be obtained, then 
@eminimum of six channels are required for effective trans- 
mission. At transmission rates very close to capacity how- 
ever, very elaborate coding scnemes are required to achieve 
@ reliable transmission. In general, it can be shown that 
mor proper coding the probability of error is bounded by 
functions of the transmission rate (32): 


271 (By, (R) CL O(n) = P(E) 2, 7nE(R) 


where E(R) and E(R) depend on the channel, and n is the code 
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block length. For rates less than channel capacity, both 


E,(R) and E(R) are positive. These functions will not be 


evaluated for this channel. 


For this fading channel, the rate is limited by the 
number of parallel channels. Tne signal interval cannot be 
shortened to much less tnan one sixth of a second because of 
intersymbol interference. The number of channels is limited 
by the useable frequency band. Above about 1500 cps, very 
little capacity is gained for the addition of an extra channel. 
Assuming a signal interval of one sixtn of a second, and trans- 
mitting five information bits during each interval, the ef- 
fect on the rate as a function of capacity of adding more 
channels is summarized in table 3, and figure 13. A minimum 
useable frequency of 900 cps and a frequency separation of 


20 cps Was assumed. 


Table 3% 
Channel # frequencies P(é) AC_ “C4 R/C 

L 900 & 920 BOOZ eG sors 

2 940 & 960 ~0040 Dies 

3 980 & 1000 0057 bey 

4 1020 & 1040 .0076 5.6 

5 1060 &1080 Ol Ow 5.5 28.4 is OG 
6 1100 & 1120 0124 54 BBG 89 
? 1140 & 1160 0154 ORS 5o.m SMe 
8 1180 & 1200 PO? le Sipe Wu 2 68 
9 1220 & 1240 nocoe 4g 49.1 Rou 
10 1260 & 1280 0333 4.7 Saae «56 
1 1300 & 1320 0408 465 58.3 so 
EZ 1340 & 1360 0562 4.1 62.4 48 
a> 1380 & 1400 0685 3.8 66.2 5 
14 1420 & 1440 .0918 3.3 69.5 43 
oy 1460 & 1480 Padi 2.9 72.4 41 


* Based on a power level allowing ten channels. For less 
channels P(€) is better, for more worse (See Part I). 


J 
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Figure 13. The effect of adding additional higher frequency 
channels on capacity and rate. (Transmission at 30 
bits per second.) 
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There are several interesting methods of coding and 
decoding which will be discussed after developing a com- 
munications system for the Nydroaceus tleme sine mea tic 
next section, the system will be developed first using simple 
transmission and tnen diversity and block coding. Block 
coding is taking a block of input digits and coding them into 
mepeonser group of coded digits for transmiss 1On vee (oserre: 
will be calculated for a decoder which makes an intermediate 
decision on each bit. As noted this is not the mosv efficient 
use of the channel, dut it is considerably easier for decoding 


ena itor error calculations. 


4, Representative Communications Systems. A simple system 


woich can transmit the necessary 30 bits per second is com- 
posed of five parallel binary channels. Each one sixth of a 
second, a five bit symbol is transmitted, one bit per channel. 
To avoid shipping noise a minimum frequency of 900 cps wil] 
be used. Table 4 lists the individual crossover probabili- 
ties for these five channels. The orthogonal signal pair are 
12 cps apart to allow a 1/6 second signal interval. The 
channels are separated by 20 cps, to permit easier tuning, 

in case of a doppler shift which could be as high as 15 cps 


more ewO Ships at 15 knots. 


Table 4 
Channel # frequencies P(E ) 
il 900 & 9Y9le 0.0014 
Z 932 & 944 0.0019 
5 964 & 976 0.0024 
Ly 996 & 1008 0.0030 
S) 1028 & 1040 ©0036 





Be 


The error rate. per five bit symbol for this transmission 
system can easily be computed as: 
P(E) = =P (now ert errors.) 

For independent bit errors the error probability is easily 
eomputed to be 0.013 per five bit symbol. Due torrecciver 
imperfections, and variation of the channel from the model 
which will be discussed in the next section, tnis can only 
be considered an order of magnitude. ‘therefore, methods of 
improving the reliability will be investigated. Shannon'*s 
coding theorem discussed in the last section indicates that. 


the error rate may be improved without an increase of vower. 


The easiest method to implement for this channel is 
Mouwenormally thought of as coding, though it is in reality 
@simple code. This is frequency diversity. For illustration 
fifteen parallel channels will be used, with each bit being 
transmitted over three channels instead of one. ‘to maintain 
the same transmitter power, each channel will use only one 
third of the power as for the simple system binary channels. 
The error probabilities for the fifteen channels are listed 
in table 5. The frequencies are spaced as for the simple 


five channel system. 


ia DiLemS 
Channel # P(E) Channel # P(E) 

ug 0025 9 0264 
2 0059 10 0324 
3 ROOT 11 . 0385 
by 0090 Le OLLG 
D 5 onlaky, 13 0550 
6 OLY“ 14 .0705 
“ POL GH 1S 0840 
8 Roee 
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A block diagram of the triple diversity receiver is snown in 
figure 14. The w; are weighting functions which depend on 
the signal~noise ratio for the particular binary channel (33): 


Ws = By /Ni 


For the fifteen channels being used, they range from 0.99 
for channel number one to 0.86 for channel number fifteen. 
For a simple equal energy L-fold diversity system the error 
probability can be computed from the following formula (34): 


oe pb S L+j-l1 J 
He HY?) a-2 


unas) ls plotted for L equal to l, 2, and 3 in figure 15. As 
an example, if p is 0.04, three-fold diversity yields a bit 
error probability of 0.00065. For a five bit symbol this 
would become an error rate of 0.0033. Assuming that p equal 
to O4 is a good average from table 5, this is a four times 
improvement over simple five cnannel transmission. A general 
error calculation method for unequal energy and noise on the 


different diversity channels is worked out in Appendix D. 


If data is to be transmitted with an accuracy better 
Maen vnis diversity system, coding will be required. One 
method, block to block coding will be discussed now, Otner 
methods will be compared to it in the next section. A prac- 
tical block to block coding system is the parity check code. 
For illustration, each five bit symbol will be encoded into a 
fifteen bit code word. The first five bits of the code word 


are the five information bits; the remainder are modulo two 
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Figure 14. Three-fold diversity receiver. 
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Figure 15. Binary receiver error probabilities Cae 
the Rayleigh fading channel. 
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sums Of several of these five bits. Froper selection of the 
code, and use of an appropriate receiver, will permit cor- 
rection of any tnree bit or lesserverror ei erecepcuc. 

of the fifteen bit code word (35). For this code the error 
probability thén bécomés: 


P() = 1l-p(no errors) - p(one error) - p(two errors) - o(three 
errors ) 


Carrying out this calculation for the data of table 5, the 
error probability is approximately 0.004 per five bit symbol, 
or approximately the same as for triple diversity. Further 
improvement can be accomplished by coding by use of a longer 
meeock lengsth, for example, coding two five bit symbols into 


Be cvnirty bit code word, or three symbols into a 45 bit word. 


5. Suitability of Codine Lechnicues. In the last section 
simple transmission, frequency diversity and a short block 
length parity check code were compared. There was a signi- 
Breantly smaller error using diversity or the coding scheme. 
To achieve a still smaller error probability, much larger 
block length or other types of codes would be required. 
Various coding techniques will now be discussed, Which may 
pe used if the added complexity can be justified by the im- 


proved performance. 


A block to, block parity .eneemeodet as used. as an 
example in the last section. For this type of coding systen, 
the number of computations required by the decoder increases 


roughly as the cube of the block length (36), and as noted, 
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Only specific error combinavions cod be corrected. IBe ois: 
limits its capability of correctine a Jarcew op -timo 1 errorc. 
except for very long codes. In general, block code receivers 
make a decision on each bit before decoding. This means that 
the confidence level in each bit is lost. The performance 
ean be improved by finer quantization in the receiver sce! ore 
decoding. It can be shown, however, that three level quanti- 
@aei0Nn Will achieve half of the possible improvement, and 
that finer quantization than tnis gives an increasing small 
improvement (37). Three level quantization is best described 
by the binary erasure channel shown in figure 16. For this 
channel the bias level must be properly set so that both p and 
q are small, where p is the crossover probability, and gq is 
the erasure probability. this type of channel is more dif- 
ficult to instrument than the simple binary channel and has 


not been used as extensively (38). 


Some of the drawbacks of block to bilock coding may 
be met with convolutional codes (39). A ie encoder is 
shown in figure 17 (40). Each transmission interval, v 
parity checks on a k bit register are transmitted. Between 
intervals, one bit is moved out of the k bit register on the 
right, and one new information bit is moved in at the left. 
@acn information bit is in the transmission sequence for k 
intervals, and tnere are an average of v transmitted bits per 
information bit. Since at least five information bits must 


be transmitted on fifteen channels six times per second, v is 


22 





Figure 16. The binary erasure channel. 
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k bit shift register 


module two adders 
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Figure Lien COonvolutionadiewcode rn: 
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limited to three for the channel being described. 


Two good decoding schemes are known for convolutional 
codes: sequential decoding (41), and threshold decoding (42). 
In the former ann bit message is treated as a tree with n 
sequential nodes. A very good reliability can be achieved 
if a computer with a large enough memory to backtrack through 
many paths along the tree is used. This may be 210 bits or 
more (39). Once an error is made because backtracking is no 
longer possible, a large block of errors may be made. This 
is satisfactory ina feedback system where repeats may be 
requested. ‘Since the transmission time over a fifty mile 
hydroacoustic channel is approximately a minute, this type of 
Channel is best suited for broadcast or prearranged time 
Sharing, with little or no possible feedback. A simpler 
decoding scheme for convolutional codes which avoids these 
block errors at the expense of a higher average error rate 


tamunresnold decoding. 


The threshold decoder uses only shift registers, modulo 
two adders, and a threshold device. For good results each 
bee must remain in the register for more than one fifteen 
bit transmission time, since for this snort coding length 
the results are about the same as for triple diversity as 
noted on figure 15 (43). For short block lengths, the results 
with threshold decoding on a convolutional code are somewhat 
meuover than for block to block coding, with no more complexity. 
This is because it can take advantage of the a posteriorl 


probabilities in the binary channels. 
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DISCUSS ION 


1. General. In this section the assumptions leading to 

the channel model and the resulting communications system 
Will be discussed first. The importance of the assumptions 
to the resulting communications system, and their expected 
Variations will be brought out. Next the hydroacoustic con- 
ditions in various ocean areas at different times of the year 


Will be discussed. 


2. Assumptions. The most important assumptions made in the 
preceding development have been the following: 

a. there is no intersymbol interference. 

be The multipath reception produces a Rayleigh 
amplitude distribution. 

ec. During each signal interval, the received phase 
and amplitude are constant. 

d. The receiver is syncenronized to the transmission 
interval. 


e. The noise is white and Gaussian. 


f. Errors on different binary channels are independent. 


Assumptions a., c., and d. are related. To minimize inter- 
symbol interference, the signalling interval must be long 
enough for most of the energy transmitted during a given 
interval to be received during an equal delayed interval at 
the receiver. The time used is based on the calculation in 


aopendix CC. The intersymbol interference can be treated as 


“a J 


J 
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an additional noise component wnich increases tne crossover 
probability, or if it can be estimated, used in a receiver 
with memory. By making the signalling interval long enough 
for intersymbol interference to be negligible, the interval 
may become too long to assume a constant phase and amplitude 
during the interval. This problem can be corrected by using 
time diversity reception during each signal interval. Since 
a small amount of intersymbol interference is unavoidable, 
and also because the transmission time between the terminals 
Will vary, it will be impossible to obtain perfect synchroni- 
zation at the receiver. Range changes may be corrected by 

a doppler circuit. Transmission path variations should be 
fovr ly small so that the time variations due to this phenomena 
Wil be small compared to the signal interval. in tnis case 
tnese variations are also correctable. In summary, these 
assumptions can be met as well as desired by proper selection 
of the signalling interval and the use of time diversity 


meeeDLion if required. 


Experimental data indicates that the Rayleigh ampli- 
tude distribution is a good assumption much of the time, but 
that it may differ (19,20). For long ranges these differences 
are most evident in the higher moments, and do not differ 
Significantly in the mean or mean square. At short ranges 
only a few paths may be of interest, and in this case the 
distribution is more favorable than a Rayleigh distribution 
(19). Figure 18 compares a typical short range amplitude 


distribution with the Rayleign distribution. 
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Figure 18. Typical amplitude distributions. 
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meeure 19. Stfect of a sloping bottom on communications 
via the bottom reflection channel in shallow water. 
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Although most ambient sea noise can be considered 
White and Gaussian for the narrow band binary channels, other 
forms of noise may be present. Rotating machinery harmonics, 
Or marine life may affect only one or two frequencies. Most 
of these noises are most severe under 1000 cps, but a source 
on a receiving ship is potentially dangerous. A careful noise 
analysis of potential receiving vehicles will permit optimum 
frequency selection with respect to the noises. Impulse 
noises may also be present. These can be minimized by clip-~ 
ping the signal at a level which depends on the average im- 


pulse noise power and the expected signal levels. 


In arriving at } system error probabilities, errors 
were considered independent of other channels. If these 
binary errors are not independent, then the resulting system 
error probabilities will be different. If for example the 
errors on two of three diversity receivers are correlated, 
the resulting bit error is only about as good as double and 
not triple diversity. No sood data is available on how closely 


correlated transmission over adjacent frequencies is. 


3. Ocean Conditions and the Channel Model. Several factors 
affect the prevailing ocean hydroacoustic conditions (40). 
The temperature, and hence the vertical profile of the velo- 
city of sound vary with both latitude and time of year, and 
in some cases the time of the day. The depth of the water 
and the bottom profile affect bottom reflection transmission, 


meavner conditions affect both noise and surface reflection 
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losses. All the calculations tiade have been for a moderate, 
state three, sea. Heavier weather conditions reduce the 


reliability considerably. 


ourface trapped channels occur only when the water 
temperature is isothermal or perhaps increasing with depth. 
this condition is somewhat unstable and can occur only when 
the surface waters are being cooled. #xperimental data 
indicates that these conditions exist most of the year in 
polar regions, during winter in mid-latitudes, and only 
rarely in tropical waters. In mid-latitudes the layer depth 
will be from 1000 feet to 1500 feet in the winter, but a 
sharp thermocline will exist near the surface in the summer. 
Since Water depth is not important to the surface channel, 
the model developed will apply to all high latitude ocean 
areas in the winter, but only to snallow Water, or polar 


Waters during the summer months. 


The bottom reflection channel can also be quite 
Variable as to results. For any given depth, the bottom 


slope is quite important, as is illustrated in figure 19. 


ey / 
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Pare sy 
CONCLUSIONS AND RECOMMENDATIONS 


1. Conclusions. In Part I it was shown that the signal- 
noise ratio for the surface channel and the shallow water 
bottom reflection channel could be estimated. It was further 
shown that the appropriate model was a Gaussian random pro- 
cess, in Wnich the received signal had a Rayleigh amplitude 
distribution and unknown phase with additive white Gaussian 
noise. Using this model an appropriate binary receiver was 
found to be a quadrature correlation device. A combination 
of five of these receivers could be used to send data at 60 
Wpm With a moderate error rate. Improvement in the error 
rate Was best obtained by frequency diversity, and could be 
further improved by coding at the expense of complexity and 
cost. In Part III, it was shown that the model was reasonable 
for shallow water any time if tne bottom did not slope away 
from the transmitter too fast, for mid-latitudes in the 
winter, and for polar regions most of the year. these 
results are applicable for moderate weather conditions, and 


Bemees Of up to fifty miles. 


2. keconmendations. Although the basic form of the receiver 
could be developed, the actual parameters could only ode 
estimated. A better understanding of these parameters is 
needed for tne actual development of a suitable systen. 


Some items of interest which warrant additional experimental 
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research are the following: 


the noise Spectrum Cf povem pia mec iy wie 
Mon Le wea. 

the range and distri out Pemeon side oo ee ae 
delays for various channels. 

the length of time during which phase and ampli- 
tude of a multipath signal will remain coherent. 
the correlation of errors on binary channels of 


Sa iacenG 1requcicy., 


ay 


\! 





Surface Channel SamplesGr boule ceqee 





these calculations the following assumptions will be made: 
Isothermal water with a velocity gradient of 0.0052 m./ 
sec./m. 

Prequence mea Lee 

fManee ... 50 miles. 

Bayer cepth ... 30.5 meters. 

sea State ... 3 

Source ... 185 db (ref. 0.00002 newt./m,) at 1/4 nm. 
medaius. inis is 173 db at a 1 metver wadius. 


a = Cs 


aS CT (sin 9 - sin 9) 


a= 0 
6 = 1500 m/sec. 
a = 0.0052 m/sec/n. 
me— 90,5 1 
cos 9, = Ce ~ TO 
c, taY Zz 
Pee OU.0152 radians = 0.07 dezrees 


R = 4380 meters 

Spherical spreading loss. 

L = 20 log R/2 = 20 Were ae al)0. = 67 db. 
Cylindrical spreading loss. 


L = 10 log 50 miles = 10 log Siow PIS tellon 
R/2 


pbrrace re: lection Loss. 


Number reflections = 50 miles = 9.2 
Passi 





ia 


Vs 
Og 


Number reflections = 10 

Loss per bounce = 0.95 db for a 4 foot wave height 

L = 0.95 X 10 = 9.5 ab. 

“absorption. Negligible for this range and frequency. 
Total losses. 67 + 15.7 + 9.5 = 92 db. 


Sienal levele. 173 = 92 = 81 db (er. 0700002 e.e 1. a 


Moccon Reflection Cranneleealcuteuro oe 


Tne same assumptions are made as for the surface channel 


except there is no layer, and: 


Peacer depth 2... 300 meters, fiav boreomn, 


WeloeouLy eracient <5. 0s0/5 me/Ssec, mu. 


H= Co (sin Q) 
acoso, 
me= 3500 meters 


Spnerical loss. 


[Le Oe ejes 


jo 


oO = 71 doe 
eylindrical loss. 


a= 10 10%» SO nae = 10 loo 230i cma 


mouLom Loss. 


Number reflections = 50 miles = 11.5 
ZR 


Number reflections = 12 

ieee 0.3 = 3.6 db. 

Suri aces Oss . 

Ween 95 = 11.4 dab. 

Rosmbosse 71 + 13.6 + GlOmmee) — 100 db. 

Signal level . 173-100 = 73 db. (ref. 0.00002 newt./n.) 
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ie Cavitation. 
Cavitation occurs when the acoustic pressure plus the ambient 
pressure is less tyan the veanor pressure of Water. Tne 
following numerical values will be used in this calculation: 
vapor pressure = 0.6 vsi at 85 degrees F. (decreases with T) 
atmospheric pressure ... 14.7 vsi 
mecver depun = 20 feet 
a. otatic pressure. 
Dee Lo7 S 0.444 X 20 = 23.5 psi 
pee Sllowable acoustic pressure. 
mee) Zou nO.O0 = 22.9 psi = 1258 x 16 SO: 
c. Sound pressure level. 
pet = 20 log p/pi(ref) 
o(ref) = 0.00002 newt./m. 
SPL = 20 log,7.9 X 10 = 198 ab. 
Allowing 3 db for peak versus mean square, the SrL is 


limited to 195 ab (ref. 0.00002 newt./m.) 


2. Non-Linearity. 


For the acoustic approximation we require: 


eq 
S 


WIU ¢& 
Ix 


q 
ct 


Assuming spnerical harmonic waves, the following equations 


ea ly: 


H 


Do wet —| 
= ej ™t-kr) 


p 


eaves ped (t-kr ) 
(2+ jie) = 





Taking the indicated derivatives we now obtain: 


5% a yee 
ie jaa oe] « ) We 


sf Comnliet = mate 
_— Ld, —_ 
bet. t, = are = OO mene 


ca =——lemecer 
° li. det — ‘ 
Baen, Since kK = =, the tracketed Cerm on) these? t ts less 


than 10. ‘the inequality can be simplified to 


2 
2: K< We 


In sea water, e, = 1000 kg/m. 
2 6 
w =2hrxXi1o atf=1 ke. 
Beecnat the inéguality becomes: 
z q 
p«<47 X 10 
8 
Letting p = 10 , We obtain the following sound pressure level: 
SPL = 20 log, p/p(ref) 


rg 2 
= 20 log 10 = 254 db (ref. 0.00002 newt./n.) 
"0.00002 
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1. time Disnersion. The maximum difference in the length 
of time it takes a signal to reach the receiver via two dif- 
ferent multipaths limits the rate at which Sienals may be sent, 
This difference is due to sound velocity fluctuations, and 
math length fluctuations. nis time difference will now be 
estimated for tne special case of a bottom reflecting channel. 
Variations will be less severe for surface channels or for 
shallower Water. It will ce assumed that tne averazse velocity 
of sound differs by a maximum of four meters per second over 
different paths. This corresponds to an average temperature 
difference of 1.0°F. This is a reasonable value since va- 
riations over several miles are normally quite small in open 
waters. It will furtner te assumed that the ray patns are 
straight lines. ‘this is illustrated in figure 20. Since 9 
femlimited to 10° - 15° at the bottom for sood reflection 
coefficients, an average value of 6° will be chosen. This 

is reasonable since the most significant paths are those 
leaving the transducer within a few degrees of the horizontal. 
All these paths will reach the bottom at avvroximately the 


same angle, so that the average Variation Will be onlyeeccurel-s 


Using these values, a ranze of 50 miles, and a sound 
velocity of 1500 m/sec., the time variation will now be 


computed. 


ah 


D/V = slant range/velocity 


R 


DeCOswer=  1OC1ZOnlLal range 
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Higure time dispersion GCalculaciom: 
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dv - tan @de 


eile 
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<j lr 


2 
< [0 
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R= 50 miles = 80,500 meters 
¥ = 1500 m/sec. 

dV = 4m/sec. 

© = ,104 radians 


eo = 0174 radians 


Carrying out the calculation the time variation is 0.240 sec. 
7nis is a Variation from tne straight line path, however. 
maven the shortest patns will range from the straignt line by 
-0OF - .06 sec. This leaves a total variation of about 0.2 
seconds which compares favorably with the only experimental 


Meee available (23). 
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ArrENDIX D 


i. Diversity Error Calculation oo ce) eo iG Lee 





tiat the diversity receiver makes a decision that m, is 


Sent if: 
2L 2L : 
a vy, Qa 
Wie “Le z = Ai Al 
42} iia 


macre the ay; are statistically independent zero mean Gaussian 
random variables. Since the wy, are all approximately equal 
they will be neglected and the error probability will be 


computed for making the decision on: 
ra 
Re a. > g a: 
. 4 ae 
it! esl 
Ey pairwise addition of tne quadrature components on each 


Channel, the decision can also be represented as: 
L 
a 


L 
gh > be 
ipl : 


4 
= 


4 


where the b; and the b, are svatistically independent Rayleigh 
eietributed random variables, with tne following probability 


eso curicutions. 


2 





: = Lin. Z i < ._< 
P(; a aa Sore, aS 5, 
=a) ted) cee 
p(b; ) = ae Roe 2b: 0.4 bass 
‘ N; 


é 


The error probability can be easily calculated if the vro- 
babilities of tne two sums are known. Tnese will now be 
out ina closed form. Define a new variable: 

2 
then it is easily shown that: 


p(k) ) = _ 27 exp - 2X: 
B, +N; BE. tN: 
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tne probability distributionvweomsseccuesumson=cac kj cen be 

found by noting that tre Ciaescecruseurc. nunciler) oOfet ewe ua 
of independent random variables is the product of their in- 
@ividual characteristic Tuncti ener eee erase teristiceou4— 


@o2L0n of k is: 


[oy ea = 1 
vd B. ae s+2/E; +N; 


' 


Z| 
d 


nerefore the characteristic function of the sum of the Kj sue 
p(aiy= 8) eS8i an, = (2) (_i_ 

° at Ea thy s+2/E, +N; ’ 

By a partial fraction expansion it is easily shown that the 


probability density function of tne sum of the xk; is the 


sum of L exponentials. 


Knowing the probability functions, the probability of 


error can easily be calculated. If m, is sent, and 


z 
so 
ea. 7 


NM 
a 
i 


R, 
meen tne error probability is: 
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This integral can be solved in general. For the special case 
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Figure 21. Diversity error with unequal signal-noise ratios. 
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